Competition between reduced delocalization and charge transfer effects for a two-band Hubbard 

model 



O 
00 



I 

O 

o 



> 

m 



o 



X 



Reza Nourafkan and Frank Marsiglio 
Department of Physics, University of Alberta, Edmonton, Alberta, Canada T6G 2G7 

We use the embedding approach for a dynamical mean-field method to investigate the electronic properties 
of a semi-infinite two band Hubbard model at half- and quarter-filling. Two effects determine the degree of 
correlation at the surface: first, there will charge transfer between the surface and the bulk, and, secondly, 
electrons at the surface are less delocalized due to the reduced coordination number We determine the result 
of these two effects and compute the quasiparticle weight. It is shown that depletion of charge from the surface 
to the bulk at quarter-filling competes with enhanced correlation effects; the net result is that at quarter-filling 
the quasi particle weight at the surface is approximately equal to the bulk quasi particle weight. Only when 
the charge transfer approaches zero at large interaction strengths does the quasi particle weight at the surface 
become lower than that in the bulk. 

PACS numbers: 71.38.-k, 71.30.+h, 73.20.-r, 71.38.Ht 



I. INTRODUCTION 

Theoretical research on strongly correlated electrons in in- 
homogeneous systems within the framework of a half-filled 
semi-infinite single-band model has shown that both the re- 
duced coordination number and (probably) reduced hopping 
elements near the surface cause an enhancement in correlation 
effects at the surface relative to the bulk. '"^ This in turn leads 
to an exponential decay of the quasi-particle weight z,^ as one 
approaches the surface from the bulk. A lower value of z at the 
surface means the coherent peak near the Fermi energy has a 
lower residue at the surface and the incoherent Hubbard band 
is more pronounced than in the bulk, in agreement with sur- 
face sensitive photoemission experiments on transition-metal 
oxides such as Cai_:cSrj;V03 andLai_2;Ca^V03.'^"" 

Recently, we emphasized the role played by charge transfer 
between the bulk and the surface, when doping a Mott insu- 
lator, and concluded that in approaching half-filling, and for 
large Coulomb interaction, charge accumulates at the surface 
from the bulk.'^ A larger charge density at the surface en- 
hances correlation effects and causes a reduced quasi-particle 
weight at the surface. Therefore, at the surface of a strongly 
correlated electron system described by a single-band Hub- 
bard model, two mechanisms lead to enhanced correlation ef- 
fects: weaker delocalization at the surface, and charge trans- 
fer between the bulk and surface layers. We note that charge 
transfer has also played an important role in understanding the 
conducting properties of interfaces between a band and a Mott 
insulator,'^ '^ and the reader is referred to the references for 
further information. 

Although in a single band Hubbard model these two mech- 
anisms both contribute to increasing correlation effects on the 
surface, their effects in a two band Hubbard model at quar- 
ter filling actually compete with one another. Indeed, due 
to depletion of charge from the surface to the bulk at quar- 
ter filling, the surface layer will have a reduced correlation, 
while reduced delocalization at the surface leads to an en- 
hancement of correlation. Therefore, the two band Hubbard 
model at quarter-filling results in a situation where the inter- 
play between these two mechanisms can be investigated. The 



purpose of this paper is to study this interplay to see if either 
mechanism dominates under various circumstances. 

The outline of this paper is as follows. In Sec. II we in- 
troduce the model Hamiltonian, which is a semi-infinite two- 
band Hubbard model with layer dependent parameters. After 
a very brief description of the dynamical mean field theory 
procedure used in combination with the embedding technique, 
we present and analyze results in Sec. III. A summary is pro- 
vided in Sec. IV. 



II. THE MODEL AND METHOD 



Our aim is to model specifically the phenomena occurring 
at the surface of a strongly correlated electron system that is 
governed by the two-band Hubbard model. The system is 
a three-dimensional, bipartite simple-cubic (sc) lattice with 
nearest-neighbour hopping only. The lattice is cut along a 
plane perpendicular to one of the coordinate axes, e.g. the z- 
axis [sc(OOl) surface]. For purposes of calculation, the system 
is considered to be built up of two-dimensional layers parallel 
to the surface. Accordingly, the position vector to a particular 
site in the semi-infinite lattice is written as RsUe — ^i + Ra- 
Here Ra stands for the coordinate origin in the layer a, and 
the layer index runs from a = 1 for the topmost surface layer 
to infinity. The vector r^ is the position vector with respect 
to a layer-dependent origin, and runs over the sites within the 
layer. Each lattice site is then labelled by indices i and a. In 
order to keep the problem simple, the interaction parameters 
are considered layer independent. To simulate different cor- 
relation strengths and relaxation processes at the surface we 
allow for a possible modification of the hopping between the 
surface and the subsurface layer ta j2 = ti2Sij. For simplic- 
ity we also assume that the hopping elements are the same for 



the two orbitals. In this notation, the Hamiltonian reads: 
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where diamcr ( di^ma I i^ ^^^ destruction (creation) operator 
for electrons with spin <t and orbital index m on site i of 
the a layer. The orbitaFs electron density on site ia is de- 
noted riiam, and tjajp is the hopping matrix element between 
two nearest-neighbour sites. U is the intra-orbital Coulomb 
interaction, U' is the inter-orbital Coulomb interaction, and 
J is the Hund coupling. The last line of Eq. (1) shows 
the pair-hopping terms. We fix the energy scale by setting 
t(^iajp) = t = 1 for a, /? 7^ 1 and we adopt the conventional 
choice of parameters, U' = U — J, which follows from sym- 
metry considerations and, for definiteness, we set J = C//4. 
The results presented here are also valid for other values of J 
(provided U' — J remains positive). 

The potential Via shows the electrical potential at the site 
ia, which comes from a redistribution of the charge near the 
surface; this potential obeys the Poisson equation. For a bulk 
lattice, in a phase with translational invariance, the local oc- 
cupation is independent of site, {niaf) + (riial) = (Uia) = U. 
On the contrary, for a semi-infinite lattice, the different local 
environment of the surface sites causes the local occupation 
near the surface to differ from the bulk filling. This charge 
redistribution on the layers near the surface gives rise to an 
electrical potential that couples to electrons and modifies the 
Hamiltonian. Due to two-dimensional translational invariance 
parallel to the surafce, the spatially varying potential is a con- 
stant for all sites in a layer parallel to the surface. Via = Va, 
and should be determined self-consistently. 

Our calculations are based on the embedding approach^'^ 
for dynamical mean field theory (DMFT)''^ for a simple cu- 
bic lattice. This method is based on partitioning of the lay- 
ered structure into a surface region which includes the first 
A^ layers and an adjacent semi-infinite bulk region (substrate) 
to which it is coupled. Then the effect of the substrate on 
the surface region is described by an energy-dependent em- 
bedding potential. The electrical potential Va in the surface 
region obeys the ID Poisson equation, whose discretized so- 
lution is 
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where n^ is the electron density in layer 7, ribuik is the bulk 
electron density and A = e^/ea with e the background di- 
electric constant and a the inter-planar lattice constant. This 



solution is obtained by iterating the discretized version of the 
Poisson equation from the bulk layers up to the surface. Here- 
after we set Viauik = 0. 

In this study, the number of surface layers is chosen to be 
N = 5 and we tested that this number provides converged 
results. Our impurity solver is exact diagonalization.^" The 
reader is referred to the references for details concerning these 
methods. 



III. RESULTS 

DMFT studies of the two-band Hubbard model for the bulk 
system have shown the existence of a Mott metal-insulator 
phase transition at any commensurate band filling, such as 
quarter-filling (n = 1) and half-filling (n = 2).^''^^ For J = 
the corresponding critical interaction Uc is maximum for half- 
filling. Including a non-zero value of J has two competing 
effects, depending on the filling: at half-filling a non-zero J 
lowers the critical interaction strength Uc, while for all other 
commensurate fillings the Uc is pushed to very high values by 
J. 

As we mentioned earlier, in a semi-infinite system the phys- 
ical quantities near surface are layer-dependent. In partic- 
ular, the surface electron occupation can differ from that of 
the bulk. We first study this phenomenon. At half-filling any 
charge modulation is excluded by particle-hole symmetry^^ 
and the local occupation on any layer, including the surface, 
coincides with the average filling, tIq = 2. Away from half- 
filling this is generally not the case. The charge redistribution 
in the surface region is driven by a narrower local density of 
states (LDOS) at the surface relative to the bulk LDOS (see 
Appendix A). However, the resulting electrostatic potential 
tends to restore the system towards a homogeneous charge dis- 
tribution. The strength of the electrostatic potential depends 
on the dielectric constant, for which an appropriate value dif- 
fers from material to material, and is often not well estab- 
lished. The two panels of Fig. 1 show the charge density 
profile in the surface region for two cases with a rather large 
charge transfer, i.e., U/t = 1, for two choices of tu/t ~ 1.0 
and ti2/t = 0.5, and for three values of X/t. The refer- 
ence case is X/t — 0, which gives the result without consid- 
ering the electrostatic potential provided by Poisson's equa- 
tion, or equivalently, the case with infinite e. The two curves 
with X/t = 0.24 and X/t = 0.9 correspond to e = 15 and 
e = 4, respectively. The dielectric constant of strongly cor- 
related metals is generally taken to be much higher than both 
these values.^"* As it is seen from the solid line in both panels, 
charge transfer diminishes with increasing distance from the 
surface and for the third layer, the charge density is essentially 
the same as the bulk density. Thus, these density changes are 
very local and we expect that they will be largely unaffected 
by considering the long-range Coulomb interaction. This in- 
tuition is confirmed by the results given by the dashed and 
dotted lines in both panels. Including the electrostatic poten- 
tial in the calculation with two typical relatively small values 
for e causes only a small suppression of charge transfer. For 
this reason we will not consider the electrostatic potential in 
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FIG. 1. Variation in layer occupation of layers in surface region for 
U jD = 1/6 and three values of background dielectric constant, a — 
1 shows the topmost layer. Top panel: ii2/i ~ 1.0. Bottom panel: 
ti2/t = 0.5. 



the results that follow. 

The top panel of Fig. 2 shows the calculated quasiparti- 
cle weight Za=i = (1 — 9S]Q.=i(a;)/9a;|„=o) of the semi- 
infinite two-band Hubbard model at T = as a function of 
interaction strength, U/D, where D is half the band width. 
Here, Sa(w) is the self-energy for layer a. The quasiparticle 
weight is a measure of the metallic nature of a system, with 
z = 1 for a non-interacting metal and z = for a corre- 
lated insulator. As expected, both the bulk and surface quasi- 
particle weights decrease monotonically as a function of the 
interaction strength, and they eventually vanish for values of 
U beyond the critical value, Uc- For any value of the inter- 
action strength, the quasiparticle weight of the surface layer 
Zsuii is significantly reduced compared to Zbuik, which can be 
understood as the effect of the reduced delocalizaton effect of 
kinetic energy and therefore enhanced effective correlations. 
The differences between the Za and the bulk z diminish with 
increasing distance from the surface and for the third layer, 
the quasiparticle weight is almost indistinguishable from the 
bulk z on the scale used (not shown). 

As is seen in the bottom panel of Fig. 2, the trend at quarter- 
filling is different from that at half-filling. In order to un- 
derstand quasi-particle weight behaviour at quarter-filling we 
should bring into account charge transfer effects. Indeed, 
charge depletes from the surface to the bulk in a range of den- 
sities including quarter-filling. In other words, while the elec- 
tron density in the bulk is ribuik = 1-0, the surface electron 
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FIG. 2. Bulk and surface quasiparticle weights z for the two band 
Hubbard model as a function of the interaction strength U/D. Top 
panel: half-filling, bottom panel: quarter-filling. 



density is a little less than quarter-filled, Uguii ~ 1.0, which 
results in reduced correlation effects at the surface. In this 
case, the characteristics of the surface quasi-particle are ap- 
proximately the same as those of the bulk quasi-particle, up 
to a relatively large interaction strength; eventually Zsmi falls 
below Zbuik- 

A summary of the situation at quarter filling is as follows. 
Reduced coordination number at the surface leads to two 
mechanisms: first, an enhanced ratio between the interaction 
potential and the kinetic energy, and, secondly, charge trans- 
fer between the surface and the bulk. The first effect increases 
correlation effects, while the second decreases them. Fig. 3 
shows charge transfer as a function of interaction strength. As 
is evident from this figure, at large interaction strength the sur- 
face charge transfer approaches zero and, for these interaction 
strengths, only the first mechanism remains. As expected for 
these interaction strengths, enhanced correlation effects lead 
to a smaller quasi-particle weight at the surface (see the bot- 
tom panel of Fig. 2). 

As a function of filling, the variation of quasiparticle weight 
is shown in Fig. 4; at quarter-filling the surface and bulk values 
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FIG. 3. The charge transfer Aua^i — Ua^i —n between the surface 
layer (a = 1) and the bulk as a function of the interaction strength 
U/D for uniform t and ti2/i = 0.5. 
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FIG. 4. Bulk and surface quasiparticle weights z for the two band 
Hubbard model as a function of the filling (between quarter and half 
filling) for two representative interaction strengths U/D. 



are relatively large and essentially equal to one another. They 
tend to spread apart as half-filling is approached, although, 
in cases where U/D is sufficiently large that the material is 
an insulator, then both surface and bulk values approach one 
another with value zero. Curves for two representative values 
of U/D are shown, illustrating this behavior. 

The presence of a surface is likely to alter the effective 
parameters entering Eq. (1), so we also investigate a case 
where the hopping between sites on the surface and neigh- 
bouring sites immediately below the surface is reduced, i.e. 
^12/^ — 0.5. This choice gives us the opportunity to study 
the effect of relaxation of the interlayer distance and also 
changes the competition between the two mechanisms men- 
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FIG. 5. Bulk and surface quasiparticle weights 2 for the two band 
Hubbard model as a function of the interaction strength U/D. Top 
panel: half-filling, bottom panel: quarter-filling. 



tioned above. For trz/t = 0.5 the ratio between the inter- 
action potential and the kinetic energy is further increased. 
One might expect that the kinetic energy is reduced, but that 
perhaps the charge transfer may be unaltered, resulting in a 
change of the balance visible in the bottom panel of Fig. 2. 
However, as is apparent in Fig. 3, a reduced hopping element 
between the surface and subsurface layers leads to a narrower 
LDOS at the surface and a larger magnitude of charge trans- 
fer. Fig. 5 shows the corresponding calculated quasiparticle 
weight for the half- and quarter-filled case. While at half- 
filling the quasiparticle weight is reduced further than in the 
uniform case (Fig. 2), at quarter-filling the two mechanisms 
continue to balance one another for weak and intermediate 
coupling strengths; only at very large interaction strengths, 
where the charge transfer is essentially suppressed (Fig. 3) 
can we see the effect of further enhanced correlation effects 
at the surface. For these interaction strengths the Zgurf ap- 
proaches zero faster than the uniform case. Actually, at small 
to intermediate interaction strengths the surface quasiparticle 
weight is slightly larger than the bulk value of z, though it is 
not evident with the scale used here. 



Note that we have confined ourselves to a bipartite lattice, 
in which there is no charge transfer at half-filling. For a non- 
bipartite lattice, it is possible to have charge transfer even at 
half-filling. Therefore, for such lattice structures we expect to 
have competition between the enhanced ratio of the interac- 
tion potential and the kinetic energy, and charge transfer, even 
at half-filling. 



IV. CONCLUDING REMARKS 

We have investigated the Mott metal-insulator transition at 
a solid-vacuum interface at zero temperature in the framework 
of the semi-infinite two-band Hubbard model at half-filling 
and quarter-filling. Using the embedding approach to extend 
dynamical mean-field theory to inhomogeneous systems, it is 
found that at half-filling, as a result of the surface narrowing 
of the local density of states, correlation effects are more pro- 
nounced at the surface than in the bulk. However, at quarter- 
filling the surface narrowing of the local density of states leads 
to an additional effect, namely charge transfer, and the subse- 
quent change in the surface occupation from quarter-filling in 
fact reduces correlation effects at the surface. It is shown that 
the interplay of these two mechanisms gives rise to approxi- 
mately the same quasiparticle weights at the surface and in the 
bulk. At very large interaction strengths the charge transfer 
approaches zero and the enhanced correlation causes a reduc- 
tion of the surface quasiparticle weight. 



in which p\ (E) is the non-interacting LDOS at site i. The 
variance in the LDOS is given by A'^pi = M^ ' — 
(A/f '°))2 = J2j_^, tiy For a site in the surface Y.j^i i% = 
4<2 + tl^ while for a site in the bulk it is Y^j^i tfj = 6*^- The 
reduced coordination number and probably reduced hopping 
elements of a site on the surface layer thus implies a reduced 
width, A^p- , for the surface LDOS. This key effect of a nar- 
rowing of the surface LDOS holds for an interacting system 
as well. 

The problem of how the narrowing of the LDOS leads to 
charge transfer is straightforward in a one-band model. ^^'^'' 
We will explain first for a non-interacting case and then gen- 
eralize results to the interacting one. If a completely homoge- 
neous charge distribution is assumed for a moment, the band 
narrowing implies different Fermi energies for the bulk and 
the surface. For n < 1 the surface Fermi energy lies above 
that of the bulk and for n > 1 we have the opposite situation 
and the surface Fermi energy lies below that of the bulk. To 
restore thermodynamic equilibrium and thus merge the Fermi 
energies, one has to allow for a charge transfer Consequently, 
band narrowing leads to na=i < ribuik for rn,uik < 1 and 
na=i > Ubuik for ntuik > 1- 

In the strongly interacting case (where an upper and lower 
Hubbard band appear), the < n < 2 domain of the non- 
interacting case maps onto the < ?i < 1 domain. 



Appendix A: Charge Transfer Due to a Surface 

The band narrowing that occurs at a surface can be un- 
derstood by referring to moments of the non-interacting local 
density of states (LDOS): 
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